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Abstract. Consider a family of planar systems depending on two parameters 
(n, b) and having at most one limit cycle. Assume that the limit cycle disappears 
at some homoclinic (or heteroclinic) connection when $(n, b) = 0. We present 
a method that allows to obtain a sequence of explicit algebraic lower and upper 
bounds for the bifurcation set $(n, b) — 0. The method is applied to two qua- 
dratic families, one of them is the well-known Bogdanov-Takens system. One of 
the results that we obtain for this system is the bifurcation curve for small val- 
ues of n sriven bv b ~ ^n 1 / 2 + -^-n 30024 3/2 2352961656 2 , n( n 5/2\ 

ues OI 11, given oy — ? « + 2401' 1 45294865 " 11108339166925 "- + u \ n )■ 

We obtain the new three terms from purely algebraic calculations, without eva- 
luating Melnikov functions. 



1. Introduction 

Consider a smooth family of planar differential equations (x, y) = (P(x, y; n, b), 
Q(x,y;n,b)), (n,b) £ M 2 , for which the existence of at most one limit cycle is 
already known and moreover all the bifurcations occurring in the family are well 
understood. For this family we could say that the Qualitative theory of ordinary 
differential equations has achieved all its goals and all the job is done. Nevertheless, 
from an analytic and quantitative point of view there remains a crucial question: 
to determine the bifurcation curves in the bifurcation diagram of the family. 

Some of these bifurcation curves are not difficult to find. This is the case for 
instance of the curves that control the changes in the behavior of the flow near 
the critical points. These bifurcation curves correspond to local phenomena and 
so, in principle, are easier to study. A paradigmatic example of this kind is the 
bifurcation curve associated to the birth of a limit cycle due to an Andronov- 
Hopf bifurcation. The curves governing global phenomena are in general much 
more difficult to determine. This is, for instance, the case of the appearance of 
homoclinic or heteroclinic connections. 

We propose a method that, in some cases, allows to obtain explicit lower and 
upper bounds in the parameter plane for the location of the bifurcation curves 
associated to these global phenomena. As illustration, we will apply it to two 2- 
parametric families of quadratic systems having at most one limit cycle. Since our 
examples are polynomial, we will consider their phase portraits on the Poincare 
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sphere, see for instance [111 Chap. 5] or [HI Chap. 3.10], and the homoclinic or 
heteroclinic connections will also take into account the critical points at infinity. 

We start by giving a description of the method for finding the relation between 
the parameters (n, b) for which this homoclinic or heteroclinic connection exists. 
As a first step we consider an algebraic curve C = {(x,y) G M 2 : C(x,y) = 0}, 
as general as possible, but having some of the geometrical properties of the phase 
portrait of the system. For instance, it has to pass through the two critical points 
of the vector field that are connected by the heteroclinic orbit; if the critical point 
is a saddle it has to be tangent to one of its separatrices; if the critical point is at 
infinity, it has to have a branch going to infinity in the direction corresponding to 
this point, etc. 

The second step consists in using all the free parameters in C to impose that, near 
the critical points, the curve be as close as possible to the separatrix connecting 
orbit. In most cases this can be done by finding a local analytical expression of this 
separatrix in the neighborhood of each critical point. There are several possibilities 
to impose that C be as close as possible to the separatrix. If the separatrix is a 
heteroclinic orbit we can fix all the free parameters of C imposing this condition 
only in the neighborhood of one of the critical points or in the neigbourhood of 
both critical points at the same time. If the separatrix is a homoclinic orbit we 
can impose the condition only on one of the branches of the orbit that arrives to 
the critical point or on the two at the same time. In the two families of quadratic 
systems studied in the next sections all these possibilities will be explored. 

Once all the parameters of the curve C are fixed, we impose one more level 
of closeness between both objects (the connecting orbit and the algebraic curve) 
by forcing a relation among the parameters (n, b) 6 I 2 . This procedure gives a 
curve B = {(n, b) e R 2 : $(n, b) = 0} in the bifurcation space. In fact, up to this 
point the method proposed is essentially the one followed in [2] to obtain analytic 
approximations to separatrices for some two-point boundary value problems. 

The last step is our main contribution and is the one that leads to explicit 
algebraic upper and lower bounds of the bifurcation curve associated to the ex- 
istence of connections between the two critical points. This is the most difficult 
and computationally involved part of the method. We have to prove that there 
are parameters (n, b) near the curve B such that the corresponding curves C are 
without contact for the flow of the vector field and moreover that there are values 
for which the flow crosses the curve in one sense and values for which the crossing 
is in the opposite direction. This will be clearer in our applications. 

Notice that the aim of our work is in the spirit of what Coppel proposed in his 
well-known paper [5] : "Ideally one might hope to characterize the phase portraits 
of quadratic systems by means of algebraic inequalities on the coefficients" , taking 
into account the results of [7] where the authors proved that there are bifurcation 
curves in quadratic systems which are neither algebraic nor analytic. Since there 
is no hope to find analytic or algebraic expressions of the bifurcation curves, we 
try to sandwich them between two algebraic curves. 

We present below our main results. The first one deals with a family of quadratic 
systems already studied in [15] and the second one with the Bogdanov-Takens 
system. For both systems there is at most one (hyperbolic) limit cycle. The value 
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b = b*(n) corresponds, for each n, to the value where this limit cycle, which is 
born in an Andronov-Hopf bifurcation, disappears giving rise to a heteroclinic or 
homoclinic connection, respectively, see more details in the next sections. 



Theorem 1.1. Let b = b*(n) be the bifurcation curve corresponding to the hete- 
roclinic connection for the quadratic system 

P(x,y) = y, 
Q{x,y) = -x + by 



with b > 0, n > 0. 
(i) For n > 14, 

(ii) For n > 5, 
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Theorem 1.2. Let b = b*(n) be the bifurcation curve corresponding to the saddle 
loop homoclinic connection for the Bogdanov-Takens system 



with n > 0. Then b*{n) - 
(i) ForO <n< 1/40, 



x = y, 

y = —n + by + x 2 + xy, 
H(y/n), where H is an analytic function. Moreover, 
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The fact that b*(n) = H{^Jn), where H is an analytic function, is an easy 
consequence of some results of [13], see Corollary 13.51 Our main contributions 
are items (i) and (ii). Notice that they improve the well-known local knowledge 
of the function b = b*(n) = f -y/^ + 0(n), proved for instance in [S] by using a 
blowing-up process and scaling of the time and the computation of a Melnikov 
function. In Section 13751 we compare them with other results of Perko obtained in 



the mentioned paper [13]. Finally, in Subsection 13.61 we prove that 

&*(l/4) € (IS, IS) - (0.364545247,0.3645452486). 

We use this information to improve the results of [9\. This subsection also serves 
to show how to use our approach when the problem that we face depends only on 
one parameter. 

At the best of our knowledge, it is the first time that this type of results are 
obtained for global bifurcation problems. 
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2. The first family of quadratic systems 

Our starting point is the following theorem: 

Theorem 2.1. Consider the system ([T]), 

f x = P(x,y) = y, 

\ V = Q{x, y) = -x + by + xy - ny 2 , 

with (n,b) G M. 2 . Then there exits a function b = b*(n) such that b*(—n) = —b*(n), 
nb* (n) > for n ^ for which: 

(i) The system has exactly one hyperbolic limit cycle if either n > and < b < 
b*(n) or n < and b*(n) < b < 0. 

(ii) Otherwise the system has no limit cycles. 

See Figured for a numeric plot of the function b = b*(n). 

0,4- 
0.3- 
0,2- 
0,1- 



-0,1- 
-0,2- 
-0,3- 

2 4 6 8 10 

Figure 1. Limit cycles of system ([T]) for n > 0. The curve b = b*{n) 
is obtained numerically. 

Essentially its proof follows from [12] and the theory of rotated vector fields 
[01 IE]- See also 0] and [10]. The property b*(—n) = —b*(n) is a consequence of 
the fact that the transformation (x,y,t) — > (—x,y, —t) changes the sign of b and 
n in system ([T]). 

The phase portraits of system ([1]) in the Poincare sphere appear in [1]. Indeed 
the phase portrait corresponding to b = b*(n), for n / 0, corresponds to the one 
having a heteroclinic connection of the two infinite semi-hyperbolic critical points 
which are on the directions y = and x — ny = 0, see Figure [2j From the above 
considerations we already know that we can restrict our attention to the region 
n > and b > 0. The origin is the only finite critical point of the system. Moreover 
the line y — 1 = is transversal to the flow associated to this system. A key result 
for knowing, once n is fixed, if some value b is smaller or bigger than b*(n) is the 
following lemma. 
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< b < b*(n 



b*(n) < b < n 



Figure 2. Phase portraits of system ([I]) for n > and b < n. 



Lemma 2.2. Letcj)(x,y) = q(y)x—p(y) = 0, withp(y) andq(y) polynomials, be an 
algebraic curve having a non-singular branch connecting the two critical points at 
infinity of system (TjQ) ; which correspond with the directions y = and x — ny = 0. 
Assume also that <j)(0, 0) > 0. If for some n > and b > it holds that 



dx 



dy 



> (resp. < 0) 



{<j>(x,y)=0} 



and the zeroes of M are isolated, then b < b*{n) (resp. b > b*(n)), see Figure^ 




M < M > 

Figure 3. Flow on the curve <ft(x, y) = 0. 

Proof. Consider for instance the case where M is positive or zero. The inequality 
M(x, y) > implies that the connected component of the curve {(p(x, y) = 0} 
described in the statement is crossed by the flow of the vector field (P, Q) from 
the region {0(x, y) < 0} to the region {<p(x, y) > 0}. Since for b > the origin 
is a repellor, by applying the Poincare-Bendixson Theorem on the sphere to the 
positive invariant region formed by the curve and a piece of the equator (which is 
invariant by the flow), we can prove the existence of at least one periodic orbit of 
system (CQ). By using Theorem 12.11 we get that b < b*(n), as we wanted to prove. 
The other case follows similarly. □ 

The main goal of this section is to get good lower and upper bounds for the 
function b = b*{n). 

2.1. A first approach to the bifurcation curve. Our first result takes advan- 
tage of a direct application of Lemma 12.21 by taking the curves <f)(x, y) = as 
hyperbolas. As can be seen in Figure HI the curve b = b*{n) is quite well delimited 
by this simple approach. 
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Figure 4. Observe how the bifurcation curve b = b*(n) of sys- 
tem ([T]), obtained numerically, lies between the two dashed curves 
given in Proposition 12.31 

Proposition 2.3. Let b = b*{n) be the bifurcation curve for the quadratic sys- 
tem (fTT) , with n > and b > 0. Then, for all n > 0: 



n 



n. 



An 2 + 1 



< b*(n) < 



16n 3 



when n < 1/2, 
when n > 1/2. 



Proof. Consider (j)(x, y) = a + a\X + a 2 y + xy — ny 2 . Clearly, inside this family of 
hyperbolas there are many of them having the geometrical restrictions described 
in Lemma 12.21 By computing M(x,y), since on <f)(x,y) = 0, we can isolate x 
in terms of y, we get that M depends only on y. Moreover its numerator is a 
polynomial of degree 4 in y and its denominator is a polynomial of degree 2 which 
is a perfect square. There are several ways of forcing it to have a given sign in 
the region where the connected component of <j)(x, y) = of our interest lies. For 
instance, by taking n > 1/2, 



a 



1 + 4n 2 



a x = -1, a 2 



1 + 8n 2 
16n 3 



and b 



lQn 3 



we force this polynomial to be (y — 1) times a polynomial of degree two which is 
a perfect square. Concretely we get that 

-n *■ { / \ N (An 2 + l)(An 2 y - l) 2 

in the region {y < 1} and that the hyperbola connects the two infinite critical 
points, as shown in Figure By using Lemma 12.21 we get that for n > 1/2, 
(8n 2 — l)/(16n 3 ) > b*(n), as we wanted to prove. 
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By imposing that the polynomial in the numerator of M(x(y),y) be the product 
of y 2 and another polynomial of degree two, which is also a perfect square, we get 

(2n 2 + l) 2 (8n 2 + l) 2 (2n 2 + l)(8n 2 + 1) 

a — o Ta o , 7Ta > a l 



0-2 



2n(4n 2 + l) 4 ' 1 (4n 2 + l) 2 

(2n 2 + l)(8n 2 + 1) n 



2n(An 2 + l) 2 ' 4n 2 + l' 

Moreover M(x(y),y) = y 2 R\(y)/ Sf(y) > 0, where Ri and Si are polynomials of 
degree one in y, and Si(y) ^ in the region {y < 1}. Again by Lemma I2T21 we get 
that for all n > 0, n/(4n 2 + 1) < b*(n). 

To end the proof it only remains to show that for < n < 1/2, b*(n) < n. 
Indeed this is true for all n > 0. Notice that the straight line {y — 1} can not be 
cut by the periodic orbits of the system. Since, when b — n > 0, 

/ y —x + by + xy — ny 2 \ 6 — 2m/ + ny 2 
^ —y J = (!-»)» £ 

by applying the Bendixson-Dulac criterion we get that the system has no limit 
cycles. Thus b*{n) < n, as we wanted to prove. □ 

Remark 2.4. By imposing that the polynomial of degree four appearing in the 
above proof be a multiple ofy 4 we get another quite satisfactory lower bound for b = 
b*(n) given by the function b = (An 2 — l)/(16n 3 ). Notice that (An 2 — l)/(16n 3 ) < 
n/(An 2 + 1). 

2.2. Proof of Theorem 11.11 For simplicity, we only give all the details of the 
proof of the weaker result, which is the one given in item (i), which asserts that 
for n > 14, 

' 1 



An y 1 



< 



3n 3 



Nevertheless its proof includes all the ingredients for the proof of the other results. 
To see the main difference for proving the other two statements see the footnote of 
next page. Before starting the proof we want to comment that although it seems 
natural to improve item (iii) of the theorem going further in our computations, 
our results are on the limit of what we can do with the algebraic manipulator that 
we use (Maple) and the capacity of our computers^. 

To prove item (i) we will apply in this case the general method described in 
the previous section. Recall that when b = b*(n) there is a heteroclinic connec- 
tion between two critical points at infinity of the Poincare compactification of 
system (JTJ), see Figure [2] Recall also that the directions corresponding to these 
critical points are given by the lines R(x,y) := xQ2(x,y) — yP2(x,y) = 0, where 
P<i and Q2 are the quadratic homogeneous parts of P and Q, respectively. In our 
case R(x, y) = xy(x — ny). The behaviour near infinity of the separatrices of these 
critical points can be easily obtained. For instance the one corresponding to the 
direction x — ny = 0, from the center manifold theorem, can be written as 

x = ty(y) = ny + ip ( — 

\y 



x We have used a computer with 32Gb of RAM and 2.4GHz. 
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with ip being an analytic function at zero of the form 



l+n 2 -nbf n 2 + l , (n 2 + 2)(n 2 + 1) - , x 

i){w) = — 1 + w + ~^^w 2 + ± — J \ V + 0(w 4 ) 



Indeed we will use this expansion at much higher order. 
Now we search an algebraic curve of the forrrjl 



. pM u + my + u 2 y 2 + u 3 y 3 + u 4 y* + ny 5 r-. . 

x = f{y) ■= —r^ = - ( : : ~. — ^ = ny + > JkV 

Qv{y) {y - 1) {v o + viy + v 2 y 2 + y 6 ) 



which is as close as possible to the searched separatrix. This is done by imposing 
that ty(y) and f(y) coincide at the highest order possible at infinity. Concretely, 
the eight free constants in f(y) are used to vanish the coefficients of y~ k , k = 
0, 1, . . . , 7 of ty(y) — f(y), obtaining the asymptotic expansion 

where 

(-l-n 2 + nb)N 8 (n,b) 
D 8 {n,b) 

being 

N 8 (n, b) =4n 19 b - 23n 18 b 2 + 70n 17 b 3 - 87n 16 6 4 + 36n 15 b 5 - n 18 - 231n 17 b 

+ A75n 16 b 2 - 143n 15 6 3 - 100n 14 6 4 + 59n 16 - 2488n 15 6 + 3668n 14 6 2 
- 886n 13 6 3 - 232n 12 6 4 + 567n 14 - 11151n 13 6 + 10886n 12 6 2 - 1284n n 6 3 
+ 256172 12 - 28350n n 6 + 16332n 10 6 2 - 272n 9 6 3 + 7718n 10 - 44884n 9 6 
+ 12264n 8 6 2 + 17160n 8 - 44920n 7 6 + 3648n 6 6 2 + 27000n 6 - 26112n 5 6 
+ 27024n 4 - 6624n 3 6 + 14976n 2 + 3456 

and 

D 8 (n, b) = - 9n 9 b + 9n 8 b 2 - 33n 7 b + 8n 6 b 2 + 2n 8 + 16n 6 
- 52n 5 6 + 50n 4 - 28n 3 b + 60n 2 + 24. 

Therefore, given n = n, some of the values of b satisfying N 8 (n,b) = could be 
a good candidate to approximate b*(n). It is not difficult to prove that one of the 
connected components of N 8 (n,b) = has an asymptotic expansion at infinity of 
the form b = j- + + • • • . From this expansion, and taking into account the 
results of Proposition 12.31 we decide to fix 

1 a 

An n 6 



To prove statements (ii) and (iii), the only difference is that we take a rational function 
with the same structure but having higher degrees in both the numerator and the denominator. 
Concretely, having respectively degrees 6 and 5 for (ii), and 8 and 7 for (iii). 
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for some a to be determined afterwards and the Ui = Ui(n,a) and Vi = Vi(n,a) 
obtained above. Taking 4> n ,a(x,y) = q v (y)x — p u (y) in Lemma f2T2l we get that 

M n Af(y),y) =t(n,a) J^i=o r ^ a > 1 where 2 = i_ yi (3) 

z (Ei=o Si(n,a)«*) 

and /(y) = p u (y) / Qv(y) ■ In the above expression all the functions rj(n, a) are 
polynomials in n of degree 40, and all Si(n,a) are also polynomials in n. The 
function a) is a rational function, which is always negative, and the coefficients 
of n m for each ri(n, a), for i — 0, ... ,4 are, respectively, 

5 a 3 3 3a 1 3a 

~256 ~ 4' ~64 ~ ~128 ~ ~2~' 64 ~ 256 ~ 4 ' 

When a = 1/8 (resp. —1/8) all the expressions of the above list are negative 
(resp. positive). Hence we can assure that for this value of a there exist rig 
such that for all z > and all n > n$ (resp. n > ) the numerator of ([3]) 
is positive (resp. negative). Similarly it can be proved that the denominator 
never vanishes on the same region. Indeed it can be seen that the biggest value 
of n which vanishes some of the eighteen functions, rj(n, ±1/8), i = 0, . . . , 4, and 
Sj(n, ±1/8), i = 0, . . . , 3 corresponds to a zero of a factor of r 4 (n, 1/8). This factor 
is 

2560n 16 - 446720n 14 - 2294624n 12 + 1116256n 10 + 33662656n 8 

+ 98872176n 6 + 131940378n 4 + 85765842n 2 + 21990713, 

and its biggest zero is approximately 13.397. In short, by using for instance the 
Sturm algorithm, it can be proved that for a = 1/8 (resp. a = —1/8), n > 13.4 
and z > it holds that the expression in (jHJ) is positive (resp. negative). By using 
Lemma 12.21 the theorem follows. 

Remark 2.5. It is clear that our method to choose the curve <p(x,y) = 0, which is 
one of the key points of our approach, that consists in imposing that it coincides as 
much as possible with one of the separatrices, can be done by using the separatrix of 
the other infinite critical point. Also it can be imposed that the curve <p(x, y) — 
coincides in both extremes, i.e. simultaneously as much as possible, with both 
separatrices. Although we obtain also other bounds for b = b*(n), the ones given 
in Theorem \l.l\ are, at least for n big, better than the ones obtained by using these 
slightly different approaches. 

3. The Bogdanov-Takens system 

For the Bogdanov-Takens system all its bifurcation diagram, including the num- 
ber of limit cycles is known, see for instance [3l [U [10l [13]. We summarize these 
results about limit cycles in the next theorem: 

Theorem 3.1. Consider the system (j2J), 

f x = y, 

\ y = —n + by + x 2 + xy, 
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with n and b real numbers. Then it has a limit cycle (which is unique and hyper- 
bolic ) if and only if 

n > and b*(n) < b < y/n. 

Here the function b* : IR + — > R gives the curve b = b*(n), n > 0, where the phase 
portrait of (J2J) has a homoclinic loop, which is hyperbolic and unstable. Moreover 

5 

b*(n) = -y/n + 0(n), n ~ 0. 
See Figured for a numeric plot ofb = b*(n). 




(H — i — . — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i 

0,5 1 1,5 2 

n 



Figure 5. Limit cycles of the Bogdanov-Takens system. The curve 
b = b*(n) is obtained numerically. 

Similarly that in the previous section our main goal is to give analytic informa- 
tion on the function b = b*(n) above. Following [3, 8] we know that, fixed n > 0, 
a hyperbolic unstable limit cycle is created when b = y/n, via a Hopf bifurcation, 
that this limit cycle increases size when b decreases and disappears for some value 
b = b*(n), for which the system presents a homoclinic unstable saddle loop. This 
behavior of the limit cycle is due to the fact that the system is a semi-complete 
family of rotated vector fields with parameter b, see for instance [61 [12], see Figu- 
re El The uniqueness and hyperbolicity of the limit cycle is proved for instance 
in PH- 




b < b*(n) b = b*(n) b*(n) < b < y/n b> y/n 



Figure 6. Saddle loop and Hopf bifurcations for the Bogdanov- 
Takens system. 
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3.1. A first approach to b*(n). When n > 0, instead of working with the ex- 
pression ([2]), by using the translation x\ = x — s/n, y\ = y, after dropping the 
subindexes, we obtain the equivalent system 

x — Vi f 4) 

y = 2^/nx + (b + \/n)y + x 2 + xy. v ' 

A simple first estimation of the function b = b*(n) is given in next lemma: 

Lemma 3.2. Let b = b*(n) be the value corresponding to a homoclinic saddle loop 
for system ([2]), then 

max(— \fn, ^fn — 1) < b*(n) < y/n. 

Proof. We work with the expression (j3J) of the Bogdanov-Takens system. Recall 
that the limit cycle is born, via a Hopf bifurcation, when b = y/n and b decreases. 
Hence, by using the non-intersection property of the limit cycles of a semi-complete 
family of rotated vector fields it suffices to prove that, when either b = y/n — 1 or 
b = — \/n, system (Jl]) has no limit cycles. When b = ^Jn — 1, note that the line 
x + y = is invariant by the flow of (j3J), because 

dL( Q^ V + 9L( Q y ^ { 2 Vnx +{b + ^)y + x 2 + xy) = (2^/n~ + x) L(x, y), (5) 

where L(x, y) — x + y. In other words, L(x, y) is an invariant algebraic curve with 
cofactor K(x, y) := 2y/n + x. 

Assume that, when b = \fri— 1, system (jl]) has a limit cycle T = {(x(t),y(t))}, 
with period T. Since the divergence of the system is b + ^Jn + x = 2^fn — 1 + x, 
the characteristic exponent of T is 

k \= f (2Vn- l+x(t)) dt. 
Jo 

On the other hand notice that from (j^D we get that 

f T , r- , ^ f T &Hx(t),y(t)) 

/ (2V^ + x(t ))dt= * \) >' y ) >J =0. 
Jo Jo L{x{t),y{t)) 

Hence k = —T and so T would be a hyperbolic stable limit cycle. This fact is in 
contradiction with the known results, see for instance Figure [6j Hence the limit 
cycle does not exist in this case. 

To end the proof let us show that b*(n) > —\fn. It is well-known that the 
stability of a hyperbolic saddle loop is given by the sign of the divergence at the 
saddle point, see [U Chap. XI]. In our case, the divergence at the saddle point 
is b*(n) + y/n. For the Bogdanov Takens system (|2J), when the loop exists it is 
hyperbolic and unstable. Hence it holds that b*(ri) + sjn > 0, as we wanted to 
prove. □ 

3.2. Proof of Theorem 11.21 (i). We introduce the new parameter m = n 1 / 4 and 
we consider the equivalent system given in (jl]), 

(6) 



2m 2 x + {b + m 2 )y + x 2 + xy. 
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We will prove our result in a constructive way. Indeed our first steps could 
be skipped but we believe that they are useful to understand the way we have 
obtained the approximations for b*(n). 

System ([6]) has a focus or a node at the point (—2m 2 , 0) and a saddle point at 
the origin. The linear approximation to its separatrices at the origin is given by 
the two lines: 

C 2 {x, y) '■= — 2m 2 x 2 — (6 + m 2 )xy + y 2 = 0, 
or equivalently, y = a\x, where a\ is one of the values 



, b + m 2 ± J(b + m 2 ) 2 + 8m 2 
a\ := ^ . (7) 

As in the proof of Theorem 11.11 we start by finding the first terms in the Taylor 
expansion of the separatrices of a critical point. In this occasion the critical point 
is the origin. In its neighborhood we can write these separatrices as 

oo 

y = ^j(x) := ^a k x k , 

k=l 

where a± is one of the values af, or in other words, C^l, cti) = 0. The next terms 
can be found recurrently by imposing that 



d(y-V(x)) d(y-V(x)), , „ 2 , , N 
Ky dx K JJ y + Ky — v ;; (2m 2 x + {b + m 2 )y + x 2 + xy) 



= 0. 



We have, for instance, 



-1 - a l a 2 (2a 2 - 1) 



a 2 = 1 ^ ; n, «3 



6 — 3ai + m 2 ' b — 4ai + m 2 ' 

03 (—1 + 5a 2 ) 3a3 2 — + %a 2 a^ 

4 b — 5ai + m 2 ' 5 b — 6a\ + m 2 

We remark that b — ka\ + m 2 never vanishes for k > 1. 
Now we consider a cubic algebraic curve of the form 

C(x, y) := C 2 (x, y) + c 3j0 x 3 + c 2A x 2 y + Ci >2 £?/ 2 + c 0:3 y 3 = (8) 

and we search the four free coefficients by imposing that {C = 0} be as close 
as possible to the separatrices of the saddle point. In particular notice that the 
quadratic terms of C imply that this curve is tangent to both separatrices. This 
can be done by imposing that the function 

F(x) := C(x, = h( m > b ) xk > 

k=3 

be as flat as possible at the origin. By choosing suitable Cjj in terms of b and m 
we get that 

F(x) = f 7 (m,b)x 7 + 0(x s ). 

After some cumbersome computations, done with an algebraic manipulator, we 
obtain that all the solutions of the equation f-r{m, b) = are contained in the 
algebraic curve 
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(2 + 6) (6 + 1 - m 2 )(b 2 + 9m 2 + 26m 2 + m 4 ) 2 x 

(3726m 20 + 21951m 18 + 251106m 18 - 671220m 16 + 68 4 06 2 m 16 

+ 1095126m 16 + 1694526 2 m 14 + 874806 3 m 14 - 3005247m 14 

- 18827016m 14 - 156 3 3 06 2 m 12 - 401766 3 m 12 - 50460486m 12 
+ 17286120m 12 + 34 2 06 4 m 12 - 703 5 3 96 3 m 10 - 4327026 4 m 10 

- 2360050m 10 - 10215036 2 m 10 - 390741536m 10 - 4762 86 5 m 10 

- 8481126 3 m 8 + 318433006 2 m 8 + 640000m 8 - 53 6 46 5 m 8 

- 13291926 4 m 8 - 680406 6 m 8 + 91212406m 8 - 118 1 76 36 5 m 6 

- 29 1 606 7 m 6 - 17920006m 6 - 2387886 6 m 6 - 28 9 28 1 36 4 m 6 

- 1143 2 5 8 86 2 m 6 - 189 3 8 4 66 3 m 6 - 89 5 4 46 5 m 4 + 125 4 4 06 2 m 4 
+ 22 3 26 7 m 4 + 568626 6 m 4 + 42471406 4 m 4 + 24306 8 m 4 

+ 44297206 3 m 4 + 9748596 6 m 2 + 56 7 06 9 m 2 + 2439506 4 m 2 
+ 549996 8 m 2 + 30710 76 7 m 2 + 86 38 1 56 5 m 2 + 12966 10 
+ 136086 9 + 429846 s + 39 06 7 + 1140 06 6 ) = 0. 
The branches passing through the origin of the algebraic curve given above are 

7 2401 2401 v ; 

Hence to continue our study we decide to fix 

5 72 

6 = -m 2 + ^jjj^ 4 + " m 5 , (9) 

with a to be given afterwards. 

With 6 given in ([9]) we will try to fix a in such a way that the curve C be without 
contact in the negative half-plane {(x, y) : x < 0} where the saddle separatrices 
lie. To impose this constrain we compute the resultant of C(x,y) and 

dC(x,y) dC{x,y) , 2 2 2 , 

D(x, y) := - y H (2m x + (6 + m )y + x + xy , 

ox oy y / 

with respect to y. This a very huge task. After many computations^, done again 
with an algebraic manipulator, we obtain that 

P(x; a, m) := Res(C, D, y) = x 10 m 36 (r (a, m) + r 1 (a, m)x + r 2 (a, m)a; 2 ) , 

where for i = 0, 1, 2, it holds that 

rj(a, m) = r (a, m) + r\{a, m) A(a, m), 



3 Notice that this is a main difference between the proof of this theorem and the one of 
Theorem ll.il because there the algebraic curve considered is of the form x = f(y), and so the 
problem goes easily to a one variable problem, while in this case we have to evaluate the resultant 
of two polynomials. 
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where 

A(a,m) :=(5764801a 2 m 8 + 345744am 7 + 5184m 6 + 19765032am 5 

+ 592704m 4 + 16941456m 2 + 46118408) 1/2 

and all the functions r^(a,m) are polynomials in a and m. Solving the quadratic 
equation 

p 2 (x; a, m) := r (a, m) + ri(a, m)x + r 2 (a, m)x 2 = 0, 
we obtain the solutions 

-ri(a,m) + a/ ?" 2 (a, m) — 4r (a, m)r 2 (a, m) 



x = x (a, m) 



x = x (a, m) 



2r2(a, m) 

/ 63 16807 r 282475249 2 \ 2 _ 

+ V2 a + a 2 m 2 ± 0(m 3 ) 

V 25 18450 326786400 J K ' 

—ri(a,m) — ^rf(a, m) — 4r (a, m)r 2 (a, m) 



2r 2 (a, m) 

247107 . 
+ O(m). 



304580 

Notice that the two roots of the coefficient of m 2 in x + (a,m) are 
13284 r- 4428 



a = a 



;\/2 ~ 1.12, a = a" := — V2 ~ -2.61. 



16807 2401 
Hence, taking for instance a = ±1 1/4 = ±2.75 we obtain two values of a for which, 
for m small enough, the two roots of the quadratic polynomial in p 2 (x; a,m) are 
positive and so the curve C is without contact, for m small enough, in {(x, y) : 
x < 0}. To know until which values of m this last property holds we fix a = ±11/4 
and study the signs of the functions rj(±ll/4, m), i = 0, 1, 2. 
As an example we give some details for 

r 2 (-ll/4,m) = r°(-ll/4,m) ±r^(-ll/4,m)A(-ll/4,m). (10) 

It can be seen that the degrees of r^— 11/4, m) and r^— 11/4, m) are 180 and 184, 
respectively. Notice that the above expression does not change sign when 

(r 2 °(-ll/4,m)) 2 - (ri(-ll/4,m)) 2 (A(-ll/4,m)) 2 (11) 

does not change sign. Hence we have to study the zeroes of this function which 
is a polynomial in m. This can be done analytically by using its Sturm se- 
quence. The smallest positive root of the polynomial ffTTl) which also vanishes ffTO]) 
is m ~ 0.40289. Hence we have proved that for m e (0,4/10), it holds that 
r 2 (— 11/4, m) < 0. Similarly we can prove that, on the interval (0,4/10), 

r (ll/4,m) < 0, n(ll/4,m) > 0, r 2 (ll/4, m) < 0, 

r (-ll/4,m) < 0, ri(-ll/4,m) > 0, r 2 (-ll/4, m) < 0. 

As a consequence, the curve C with b given in 0, a = ±11/4, and m G (0,4/10) 
is without contact in {(x,y) : x < 0}. It is easy to prove that for m sufficiently 
small the curve C has a loop in the half-plane {(x, y) : x < 0}. Moreover with 
the same type of algebraic methods that we have used above it can be seen that 
the loop also exists in the interval m G (0,4/10). Furthermore, in this half-plane 
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the loop is crossed transversally by the flow of (Q, inwards when a = —11/4 and 
outwards when a = 11/4. Hence we have proved that when 

o = b m := — m H m ± — m , 

K 1 7 2401 4 

and m G (0,4/10) it holds that b~(m) < b*(m) < b + {m). From these inequalities 
the statement (i) of the theorem follows by noticing that (4/10) 4 = 16/625 > 1/40. 

3.3. Proof of Theorem ll.2L (ii). Clearly the ideas of the method used to prove 
Theorem ll.2[ (i) can be applied when the curve C given in ([8]) is taken of higher 
degree. Consider now the algebraic curve 

C 4 (x,y):=C 2 (x,y) + ^ r,,rV 0. (12) 

3<i+J<4 

By employing the same procedure applied to the cubic algebraic curve we obtain 
the following relation between b and n = m 4 : 

, 5 1/2 72 30024 3/2 2352961656 2 5/2N 

& = -n 1/2 H n n i/2 n 2 + 0(n 5/2 ). (13) 

7 2401 45294865 11108339166925 v ; K J 

By using the variables n and b, the necessary calculations to obtain analogous 
results, for a given interval [0, no], to those obtained in Theorem 11.21 (i) for the 
cubic curve are beyond our computational capacity. In particular we can not 
obtain a compact expression of the resultant that would give a proof that C4 is 
without contact in the half-plane {(x, y) : x < 0}. Nevertheless, once n is fixed 
we can perform all the computations and prove for instance that for n = 1/20, 



5 1/2 72 30024 3/2 . , 

-n 1/2 H n n A ' 2 - b*(n) 

7 2401 45294865 v ' 



< 10-V/ 4 , 



and for n = 1/100, 

< 10-V/ 4 . 



5 1/2 72 30024 3/2 2352961656 2 N 

7 2401 45294865 11108339166925 v ' 

In any case, we will prove that the expression (fT3"|) gives us, at least locally, the 
function b*(n). This can be done by introducing the new variables B > and 
M > as 

= b + m 2 m2 = (b + m 2 ) 2 + 8m 2 

This change of variables is motivated by equation ([7]). Notice that by using them, 
af = B±M. Now we repeat all the procedure developed in the previous subsection 
but with the curve C^x^y) = given in fTl2l and these new variables. By forcing 
the curve Ci(x,y) = to coincide as much as possible with both separatrices of 
the saddle point (four conditions of contact with one of them and five with the 
other one) we determine all the parameters of the quartic. After that we obtain 
the relation 

r, 3 ».2 180 », 4 2366307 . 505643614857 - ^ /71 .. 

B =-M 2 M + M M + 0(M 9 ), (15) 

7 2401 90589730 44433356667700 v ;> v ; 

that ensures one more level of closeness between one of the separatrices and the 
algebraic curve, for M sufficiently small. 
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Fortunately, with these new variables it is now possible to explicitly get the 
resultant P(x; M,B), between C 4 (x,y) and 

dC 4 {x,y) dC 4 (x,y) . 2 2 2 . 

1)4(3;, 2/) := ^ U H ^ (2m x + (0 + m )?/ + x + xy , 

ox oy ' 

with respect to We obtain that 

P(x; M, B) := Res(C 4 , D 4 , y) = x 14 (s (M, B) ^ T^M, B)x" 

\ i=0 

where So and Ti, % = 0, . . . , 5 are rational functions in S and M. To have an idea of 
their complexity we introduce the following notation: we will say that a rational 
function S(M,B) is of type {i, j}/ {k,l} if after simplifying it, its numerator has 
monomials of degree between i and j and its denominator between k and /. Then 

S is {0,1}/{114,132}, 

T is {91,120}/{0,0}, Ti is {109, 143}/{20, 24}, 

T 2 is {124,162}/{37,44}, T 3 is {137, 179}/{52, 62}, 

T 4 is {137,179}/{53,63}, T 5 is {137, 180}/{55, 65}. 

Note that six of the above functions are singular at B = M = 0. Inspired by (TIB]) 
we express B in terms of M and a as 

r, „,„ s 3 ,,9 180 2366307 ,. 505643614857 , ,„ - 

B = 8(M, a) :=-M 2 M + M 6 M 8 + aM 9 , 

1 K ' ; 7 2401 90589730 44433356667700 

(16) 

where a is an arbitrary parameter to be determined later. We obtain that 

P(x;0(M,a),M) 



x 14 



= Y, P i( a ) x ' 1 

M=0 i=0 



where now all the functions Pi, i = 0,1,. ..,5 are polynomials in a with degrees 
{2, 5}, {0, 6}, {0, 11}, {0, 15}, {0, 15}, and we use a notation similar to above. We 
remark that to have the cancelations between the respective numerators and de- 
nominators in the expression of P, that allow us to evaluate it at M — 0, it is 
necessary to take B as in the expression (JIB"]) . 
The polynomial P is 

P (a) = ka 2 (296751659628833594552482011388242366232495503625a 3 

+ 719549554938315584569470362390245200816649200a 2 

+ 532161050006783873283311272385459961077760a 

+ 112647235678813465306115059636253208576), 

for some positive integer k. It has a = as a double root and three negative 
simple roots, which are approximately —0.001166, —0.000895 and —0.000364. In 
particular 

P (-2500" 1 ) < and P (2500" 1 ) > 0. 

If we fix any of the two values a ± = ±2500^ it holds that P(x; 3(M, a ± ), M) 
does not vanish if M and \x\ 7^ are small enough. This implies that shrinking, if 
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necessary, the value of M we obtain that the oval of the algebraic curve C^x, y) = 
0, which lies in the negative half plane {x < 0} and starts and ends at the origin, 
is without contact with the flow of the system (except at (0, 0)). Notice that this 
oval is born from the cusp point existing when M = and B = ^(a^, 0) = 0. 

The fact that the resultant between C4 and D4 has different sign when a = a + 
and when a = a~ induces to think that the flow crosses the oval outwards in one 
case of plus sign and inwards in the other one. To prove this fact it suffices to 
check this property, for each case, on a single point of the oval. We have chosen 
the cutting point of the oval of C±(x, y) = and the x-axis. 

The four points of the curve C^x, y) = on the 0, which is a double 

one, and 



xi (a M) =-^M 2 + —M 4 - 14661 M 6 + 1097361567 M 8 

^\\ul,xvx) 2 xvi -t 3g2 ivi 168070 ^28988713600 

57336960516777 TiylO i ( 429 _ 1360220314860156764457 \ aril , r^l yf\1\ 
3110334966739000 U280 758276772825613000960000/ ~ t ~ w \ lvl )■• 

T (n M\ = 3266440450 
^2{U.,ivi) 1212150477 

, ( 688506579126746016577 , 868315765444642387622750 \ tut , Q( ]Uf2\ 
~^ \ 52895116040087791044 357042033270592589547 ) *~ W \ 1VI )' 

Notice that the point on the loop is (xi(a, M), 0). After many computations we 
obtain that 

£> 4 (si(a,M),0) = 3718 S° 7 6000 «M 11 

_ ( 1233580648740454800 , 644635919327316538269855546 \ ajU , n(M 1Z \ 
\ 3107227739 ' 575604427043451909103 / ' w \ lvl )• 

Hence the sign on C±{x,y) = of D4 is different for a = a + > and for 
a = a~ < 0, as we wanted to prove. Therefore, for M small enough^, the value of 
B corresponding to the existence of the homoclinic loop, say B*(M), satisfies 

(3(M, -250CT 1 ) < B*(M) < /3(M, 2500" 1 ). 

Collecting all the above results and writing expression (fl5|) in the old variables 
b and n we get the proof of Theorem 11.21 (ii). 

To get an idea of how far the approximation given in Theorem 11.21 (ii) works we 
compare it with a numerical approximation of the bifurcation curve. Concretely 
we have obtained a numerical approximation, working with precision 10~ 16 , of 
b = b*(n) for n G (0, 10~ 2 ) by using a Taylor's method for solving the differential 
equation. If we denote it by b = b^ nm (n) it holds that 



5 1/2 72 30024 3/2 2352961656 2 

— n t n — n — n — in 

7 2401 45294865 11108339166925 numV ' 



< 6 x 10 



-10 



for n G (0, 10 



-2> 



To determine an explicit value of M until which the result holds we could try a similar study 
to the one done in the previous subsection, but it would be extremely long and tedious. 
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3.4. More terms in b*(ri). By applying again the procedure described in the 
previous subsections, but starting with a sextic algebraic curve 



C 6 (x, y) := C 2 (x, y)+ c i4 x % y> = 0, 

3<i+i<6 



we arrive to a relation similar to (1151) that writes as: 



tj —3nf2 _ 180 yrl , 2366307 71^6 _ 505643614857 71^8 , 121044460222851597 jitIO 
7 2401 ~ t ~ 90589730 44433356667700 ' 21794117111940173000 

_ 75409774306549331412249 /if 12 , 344552497352535858777709804917 yr 14 
25960934362572014675870000 ~ t ~ 216470837823465107132247097100000 

_ 43306460616773431694096161799928995367 j\#16 
48013140478550999259657196328505023800000 

1 91720311301427439675156623493153846098504753619 yr\% 
' 174952692015527559148011719369634259959684998000000 

_ 990748106733217809261982123885784358373281388289053276929 jyr20 
3187507017704227098311349623537514328886463857480747900000000 

1 2174773094339151212022525670857933567647566708078598123403162358137 71^22 
11614798115838021630069501971110592530190950393238633306372059000000000 

+ 0(M 24 ). 



From it and the change of variables (TBI) we obtain the following relation between 
b and n, 

h* =^n l l 2 4- 72 — 30024 3/2 _ 2352961656 2 , 161066618396136 5/2 
7 ~ l ~ 2401 45294865 11108339166925 ' 2724264638992521625 

_ 28575844096870898712 3 , 37409973403083644863711656 7/2 
1622558397660750917241875 6764713681983284597882721784375 



n 4 



1301593321483486009213262204378664 
750205319977359363432143692632890996875 

1 750633455019308628819042126726886218707352 9/2 
1366817906371309055843841557575267655935039046875 L 

_ 2188961083333347178341822657596953981848275462851032 

12451199287907137102778709466943415347212749443284171484375 ' 

1 1289326941251660725073133052778275691207040442626311930438856 1 
"T" 22685152569996135996229496037325376035529199986794205676507927734375 

+ 0(n 6 ). (17) 



,11/2 



We believe that it is an improvement of the expression of b*(n) given in Theo- 
rem (H2J(ii), with seven new terms of the expansion of b*(n), but we have not been 
able to prove this fact due to the complexity of the necessary algebraic calculations. 
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Notice that the prime decompositions of the denominators of (fTTl) have a nice 
and regular structure: 

7 = 7, 
2401 = 7 4 , 





45294865 = 


5 ■ 


7 7 - 


11, 












11108339166925 = 


5 2 


. 7 10 


■ll 2 


13, 








2724264638992521625 = 


5 3 


. 7 13 


■ll 3 


13 2 , 








1622558397660750917241875 = 


5 4 


. 7 15 


■ll 4 


13 3 


17, 






676471368198328 


45978827217843756 = 


5 5 


. 7 18 


■ll 5 


13 4 


17 2 , 






750205319977359363 


•■■92632890996875 = 


5 5 


. 7 22 


■ll 6 


13 5 


17 3 


19, 




136681790637090558 


■■■55935039046875 = 


5 6 


. 7 25 


•ll 7 


13 6 


17 4 


19 2 


23, 


124511992879137102 


■■■43284171484375 = 


5 8 


. 7 28 


■ll 8 


13 7 


17 5 


19 3 


23 2 


226851525699359962 


■ ■ ■ 76507927734375 = 


5 9 


. 7 31 


•ll 9 


13 8 


17 6 


19 4 


23 3 



where the small irregularity with the number of fives and sevens could be produced 
by some cancelations with the respective numerators. Unfortunately no regularity 
appears in the numerators. 

This regularity in the denominators of the asymptotic expansion of b = b*(ri) 
at the origin could give some clues about a possible closed form expression of this 
function but, by the moment, we have not been able to obtain it. 

3.5. Relation between our results and Perko's approach. In [13J, the au- 
thor faces the problem of the global bifurcation diagram, on the Poincare sphere, 
of a different representation of the Bogdanov-Takens system. Concretely he con- 
siders the system 

{ x = Vi (18) 
y = x(x - 1) + niy + fi 2 xy 

and, among other results, he proves the following theorem: 

Theorem 3.3 Q13J). There exists a unique analytic function h(n 2 ) defined for 
all \x 2 £ R that satisfies h'(0) = —1/7, h{—ji 2 ) = — ^(^2), o,nd max(— 1,— //|) < 
A*2^(a*2) < for all \i 2 7^ such that 

(a) System ([TBI) has a unique, hyperbolic limit cycle if and only if \i\[L 2 < and 
< < \h(/i 2 )\; the limit cycle is stable if \i\ > 0, and unstable if Hi < 0; 

(b) For fi 2 7^ 0, system (1181) has a fine focus of multiplicity one at the origin if 
and only if /i x = 0; for /x 2 < 0, a unique, stable, limit cycle is generated in 
a supercritical Hopf bifurcation at the origin of (ITS]) at the bifurcation value 
Hi — 0, and it expands monotonically with increasing \L\ until it intersects the 
saddle at (1, 0) and forms a homoclinic loop at the bifurcation value fi\ = h(fi 2 ); 

(c) System ([TBI) has a homoclinic loop at the saddle (1, 0) if and only if fi\ = h(fi 2 ); 
the separatrix cycle is hyperbolic if and only if \x 2 7^ 0; it is stable (unstable) 
on its interior for fi 2 < (\i 2 > 0). 
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Let us relate the results of the above theorem with our results. We need the 
following simple lemma. 

Lemma 3.4. When /i 2 7^ the change of variables and time 

Lb*) / \ q t 

u = TT i 2x ~ 1) > v = IhVi T = — 

transforms system (|18[) into 

, du 

u = — = v, 
dr 

, dv \x\ /x 2 (2//i +H2) . 2 , 

v = — = — - H v 4~ it + w. 

rfr 4 2 

As a consequence of the above lemma and Theorems 13.14 ll.21 (ii) and 13.31 we 
obtain: 

Corollary 3.5. (i) The following relation holds 
b* ( ^) 

M^) = -f + -^ (19) 

1 18 , 3753 , 294120207 7 _ q , 

= — A*2 + /4 lA + OUl), (20) 

7 P 2401 P 2 45294865 P2 22216678333850 P2 VP2; v 7 

where the functions b* and h are the ones defined in Theorems \3.1\ and \3.3[ re- 
spectively. 

(ii) The function b*(n) is given by b*(n) = H(y/n), for some analytic function 

H. 

Proof. By using Lemma 13.41 we get that the relations between the variables in 
system (fTSj) and the ones of system (j2j) are 

A , ^2(2^1 + ^2) 

n = — , = . 

4 ' 2 

By using them we easily obtain the proof, because the relation b = b* (n) writes as 

/x 2 (2//i + /x 2 ) = , * (A 
2 V 4 

which immediately leads to (fT9l) . From this relation and the expression of b = b*(n) 
given in Theorem 11.21 (ii). the expansion (1201) follows. Item (ii) is a consequence 



of (fi~9~]) . Theorem 13.31 and the fact that h is an odd function. □ 

We end this section with several remarks. 

Remark 3.6. (i) Expression (j20~|) improves the local knowledge of the function h 
given in Theorem Iff. 31 The proof of [13] #ia£ h'(0) = —1/7, is equivalent to the 
classical one which appears in [S] and gives the term ^■sfn in expression ( 1TB"|) . Our 
approach goes much further and it is completely different. 

(ii) By using Theorem \l."A (i) we obtain the following global result: for |/i 2 | < 
10-1/4 ~ 0.562, 



1 18 3 V, N 

~r 2 + 24oi^ - K ^ 



11 , ,4 



- ^72" I/ 12 
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Notice that this inequality proves that for e (-KH^O), 

h(n 2 ) < —yH2- 

This result is coherent with the conjecture made in [13] that affirms that the above 
inequality holds for all /i 2 < 0. 

(Hi) The inequalities of Lemma Iff. ,21 translated to system f[T8"l) read as 

max(— 1, < fahifa) < 0. 

This information is already contained in Theorem \3.3[ but our proof is different 
to the one given in [13] . 

(iv) In [13] there is another conjecture that says that /^M/^) + 1 = 0(1/ '^2) as 
ji2 — ► —00. In the parameters of the Bogdanov-Takens system (J2J) it reads as 

b*{n) = \fn —1 + f n ~ 3 ) as n ~~ > 00 . 

(v) The case [12 — includes new systems which are not contained in the ex- 
pression of the Bogdanov- Tokens system written as in <^ . For instance the case 
A*i — ^2 = corresponds to a Hamiltonian system with a center at the origin. 

3.6. A final application. In [9J it is proved the following result: 

Theorem 3.7 Q9J). The system 

x = Vi ^21) 
y = fly — ax 2 + a 2 x — xy, with a < 0, v ' 

has a limit cycle if and only if 7 < (3/ a < 1, where 7 is a positive constant. 

This constant 7 is computed numerically in that paper as 7 ~ 0.864546. We 
remark that for the values a and (3 satisfying (3 /a = 7, the system ( 1211) has a 
homoclinic loop through the origin. We will improve the results of that paper. 
We start by proving the following lemma. 

Lemma 3.8. The constant 7 defined above is 

7 = 6*(l/4) + l/2, 

where b* (n) is the function introduced in Theorem lff.il 

Proof. By applying the change of variables 

1 x y 
u = - , v = 2' T = act, 

to the system (12TI) we obtain 

du 



u = — = v, 
dr 

. dv 1 ( 3 1 . 

V — — = 7 + I ) V + M + UV. 

dr 4 \a 2 



Hence, the correspondence between the parameters of system (121]) and the ones 
of the Bogdanov-Takens system (j2]) is 

n = 1/4 and b = (3/a- 1/2. 

From these relations the lemma follows. □ 
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As we have already explained at the end of Subsection 13.31 we have computed 
6* um (77.) at several values of n, obtaining in particular that 

7 ~ &; um (l/4) + 1/2 = 0.864545247421507 =: 7 n™. (22) 

Notice that only the first five significative digits agree with that computed in [9], 
0.864546. As we will prove in Theorem 13.101 at least eight of the significative 
digits of 7 num are correct and (7 — 7 num | < 1.2 x 10~ 9 . 

Let us call 7^ := (1/4) +1/2, where b k (n) is the function obtained by adding the 
first k terms of the asymptotic expansion of b*(ri) at the origin, see Theorem ll.2[ (ii) 
and expression (IT7|) . For instance 

b\{n) 5 ' ^ ' ' 7 " 



5 1/2 



n 



-n 1 ' 2 + 
7 2401 



n 



and so on. Although we have not proved neither the validity of these expressions 
until n = 1/4 nor the validity of 6^(1/4), k = 5, . . . , 11 we can compute the values 
6^(1/4) + 1/2 and compare with the numerical approximation of 7. We obtain 



7i 



73 



6 
7 

79 = 2076 
< 2 2401 

_ 39159987 



71 



45294865 
19207287903423 



76 



22216678333850 
9421002777077246787 
10897058555970086500 
11222200726046133491344191 
12980467181286007337935000 



75 



77 

78 = 

79 = 

710 = 

7n = 



935744176562 



108235418911 
2075476618505 



2400657023927 
7562725921574 



8747634600776 
1377872021601 



1593753508852 
5020759255426 



5807399057919 



60335616104987 



66123548550000 
38798344748153 



64252511900000 
99070483500549 



79842499000000 
10350610970071 



40373950000000 
85210508369767 



86029500000000 



0. 


S57142857143, 


l7i 


7num | 


< 7.5 x 


10" 


-3 


0. 


364639733444, 


72 


7num | 


< 9.5 x 


10" 


-5 


0. 


S64556876370, 


173 


7num | 


< 1.2 x 


10" 


-5 


0. 


S64543637658, 


174 


7num | 


< 1.7 x 


10" 


-6 


0. 


S64545485251, 


175 


7num | 


< 2.4 x 


10" 


-7 


0. 


364545210070, 


76 


7num | 


< 3.8 x 


10" 


-8 


0. 


364545253274, 


177 


7num | 


< 5.9 x 


10" 


-9 


0. 


364545246497, 


78 


7num | 


< 9.2 x 


10" 


-10 


0. 


364545247569, 


Its 


7num | 


< 1.5 x 


10" 


-10 


0. 


364545247398, 


Itio 


7num| 


< 2.4 x 


10" 


-11 


0. 


364545247425, 


711 


7num| 


< 3.9 x 


10" 


-12 



Notice that the values 7^ given above approach well to 7 num , and by using again 
the results of Theorem 13.101 also approach to 7. 

In fact in the sequel we will see how the method introduced in this paper allows 
to give a concrete interval where the actual value of 7 lies. 

Recall that our approach for obtaining information of whether a homoclinic 
connection appears passes trough the construction of the two polynomials 

C k (x,y) = C 2 (x,y)+ Yl 

3<i+j<k 

dC k (x,y) dC k (x,y) 



(23) 



D k (x,y) 



-y + 



-(2m x + (b + m )y + x + xy), 



dx dy 

with k > 3, and their resultant. From a computational point of view it is much 
simpler the case where all the numbers implied in their obtention are rational 
and then C k (x,y) and D k (x,y) are in Q[x, y]. It is easy to see that this happens 
when b and the eigenvalues of the saddle point given in ([7]) are rational numbers. 
From now one we will particularize our study to case n = 1/4, that is m 2 = 1/2, 
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although clearly our approach can be adapted to all values of n such that ■sfn is a 
rational number. In our case, the rationality conditions are reduced to 

b = i and (i + l\\ i= (i.) 2 

q \q 2) \s/ 

for some integer numbers p, q, r and s. This is equivalent to find integer solutions 
p, q and t of the quadratic diophantine equation 

V + 4pq + 17q 2 = t 2 . (24) 

It is well-known how to study this type of equations. More concretely, if (p , q , t ) G 
Z 3 is a particular solution of the diophantine equation 

Ap 2 + Bpq + Cq 2 = Dt 2 , 

where A, B, C and D are integers numbers, then 

p —p(u, v) = (Ap + Bq )u 2 + 2Cq uv - Cp v 2 , 

q =q(u, v) = -Aq u 2 + 2Ap uv + (Bp + Cq )v 2 , 

t =t(u, v) = t (Au 2 + Buv + CV), 

for any (u, v) eZ 2 , is also an integer solution of the diophantine equation, because 
it holds that 

Ap(u, v) 2 +B p{u, v) q{u, v) + C q(u, v) 2 — D t(u, v) 2 = 

(Au 2 + Buv + Cv 2 ) 2 (Ap 2 + Bp q + Cq 2 G - Dt 2 ). 

For equation (124")) it is clear that (po,qo,t ) = (1,0,2) is a particular solution. 
By applying the above procedure we obtain the new solutions 

p =p(u, v) = Au 2 — 17v 2 , 

q =q(u, v) = 4v(2u + v), 

t =t(u, v) = 2(Au 2 + Auv + 17v 2 ). 

From them we obtain the suitable candidates to perform our study. We consider, 
for (u, v) G Z 2 , with v(2u + v) ^ 0, 

4t,2 _ 17 2 

b = g(u,v)= U j-t G Q, (25) 
4v(2u + v) 

for which the eigenvalues of the saddle point given in (17j), af, are 

, 2u + v ^ , Av ^ 

af = — G Q and a x = -- — — G Q. 

Av 2u + v 

The above results will allow us to prove the following lemma. 

Lemma 3.9. Fix n = 1/4 and b G R. Then there are infinitely many sequences 
of rational numbers {bj}j^N such that linij^oo&j = b and the eigenvalues af of the 
saddle point 

b 3 + 1/2 ± y/(bj + 1/2) 2 + 4 

2 

are also rational numbers. Moreover these sequences can be explicitly obtained. 
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Proof. Notice that, for u^O, the map g(u, v) given in fl25|) can be written as 

b = g(u,v) = . 

1 ' 8(f) +4 

Then the first part of the proof follows by noticing that the graph of the map 

Aw 2 - 17 

w i — > G(w) = — , toeR (26) 

covers all the real line. To prove the second part it suffices to consider any sequence 
of rational numbers {wj}j £ fq tending to one of the preimages for G of b, say 
Then, clearly, it holds that for any j G N 



J V 3> %Wj + A 2 V 

and limj^oo = 6 as we wanted to prove. □ 
Finally, we prove: 

Theorem 3.10. Let 7 be the real number introduced in Theorem \3. 7\ Then 
7 G (!!f§i,!§if§) - (0.864545247,0.8645452486). 

Notice that the length of the interval given in the above theorem is smaller than 
1.61 x 1(T 9 . 

Proof. By using Lemma [3.81 7 = 6* (1/4) + 1/2. So we will fix n = 1/4 and we 
study b* := b* (1/4), for system ([2]). For the sake of simplicity we only give the full 
details of the range of values of b* that we obtain by applying our method with 
k — 3, that is by using a cubic curve fT23l) . The proof of the theorem follows by 
taking k = 6. 

The numerical approximation (J22J) of b* gives an orientation for the actual value 
of b*. After some trials we consider the values 

19\ 89 , „, /8\ 103 



b\:=G[ — = and Vi:—G[ ,- 

3 V 8 / 368 3 \3J 228' 

where G is the function given in fT26l) . and we prove that b* G (6|, 63). Notice that 
(89/368,103/228)^(0.242,0.452). 

We will apply the procedure introduced in Section[3]for n — 1/4 and b G (63, b^) 
for searching suitable cubic curves, both having an oval through the origin and 
being without contact with the flow of the system. Moreover we will prove that 
for one of the ovals the flow goes inwards and for the other one the flow goes 
outwards. 

Fix one of the values, say b = 89/368. Then we apply our procedure to determine 
Cs(x,y), by imposing all the coincidence conditions between the cubic and the 
separatrix associated to the eigenvalue af = 23/16. Recall that due to our choice 
of b we can ensure that all our computations will be with polynomials in Q[x,y]. 
We obtain 



Ci(x,y) 



_ (16x+23^)(16j/-23a) _ 30470974207443747849 3 _ 48084904789188461109 2 
368 44286028769220680429 88572057538441360858 U 

I 300486284549883520 2 _ 78703862917780480 3 
44286028769220680429 y 132858086307662041287^ ' 
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We omit the explicit expression of D^(x,y) G Q[x, y], which has degree 4. The 
resultant of C|(x, y) and D^(x, y), with respect to y, Res (C|, y) is a polynomial 
of the form x 10 (Ax 2 +Bx+C) G Q[x], for some huge rational numbers (numerators 
and denominators with more than 60 digits), and A > 0, B < and C > 0. Hence 
the resultant is positive in {(x, y) : x < 0}. 

If we compute the resultant of C^(x,y) and dC^(x,y)/dy with respect to y we 
obtain a polynomial of the form x 2 P 4 (x), with P±(x) G Q[x]. The roots of P 4 can 
be explicitly obtained (or located by using the Sturm sequences of P 4 , in view 
of the application of our approach for k > 3). Their approximate values are 
-24.478, -9.855, x ~ -1.454 and 32.737. 

The following facts are not difficult to prove for the cubic C^(x,y) = 0: 

• Its only multiple point is the origin, which is a double point. On it, the 
curve has two smooth branches tangent to the lines (16x + 23y)(16y — 
23a;) = 0; 

• The polynomial Cg(0,y) has the as a double root and a simple nonzero 
root; 

• For each x G (x , 0) the polynomial y i — ► C^(x, y) has exactly three simple 
real roots; 

• The polynomial y i — > C^(xq, y) has a simple root and a root of multiplicity 

two; 

• For each x G [—2, xq) the polynomial y i — > C^(x, y) has exactly one simple 
real root. 

The tool that we use to prove the first item is the computation of the resultants 
between C|(x, y), dC%(x, y)/dx and <9C|(x, y)/dy. To prove the other ones we com- 
pute the Sturm sequences of the polynomials C|(x, y) considered as polynomials 
in y, with coefficients in the rational functions with numerators and denominators 
in Q[x]. This can be done, except at finitely many points, given by the zeros of 
some polynomial of x that appear during the process. In particular, we can see 
that studying the four polynomials C^(x, y) for x G {—2, x , —1, 0} we have all the 
information for y i — > C^(x, y) and x G [—2, 0]. 

The above list of properties prove that the curve C^(x,y) = has a loop which 
starts and ends at the origin being tangent to the lines (16x + 23y)(16y — 23x) = 
and is contained in the strip {(x,y) : xq < x < 0}. Since we have proved that the 
resultant Res(Cg, D%,y) does not vanish on the left hand plane {(x,y) : x < 0} 
we know that the oval is without contact for the flow of the system. By studying 
the sign of on the loop we prove that the flow crosses it inwards. Hence we 
have proved that b* > b e 3 = 89/368. 

Working similarly, for b = 103/228 and ai = 19/12 we obtain the cubic 

nul \ (12a;+19^)(12y-19a;) _ 42394240475656582327 3 _ 24891674002318104595 2 
^Zy^iy) 228 66085291294166321043 44056860862777547362 X » 

I 320066250082464600 2 _ 37723016690931312 3 
22028430431388773681 » 22028430431388773681 » 

and we prove that b* < b% = 103/228. 

Before considering greater values of k we want to comment that, also keeping 
k = 3, but taking rational numbers with big numerators and denominators, we 
could improve a little bit the knowledge of the interval where b* lies. For instance 
we can prove that 
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b* e (G (§§) , G = i|) ^ (0.2525, 0.4305). 

By applying our method for k > 3 we obtain that 

b * e (li§> - (0-364312, 0.365103) for k = 4, 

& * e (Hfii, Si) ^ (0-3645448, 0.3645454) for fc = 5, 

& * e (S. Hii) - (0.364545247, 0.3645452486) for fc = 6. 

We notice that the theorem follows from the result when k = 6. 

All the computations for the six values of b to be studied follow a similar proce- 
dure to the one described above. The main computational difficulties appear for 
proving that the algebraic curve of degree k has a transversal oval passing through 
the origin due to the high degree and huge size of the coefficients of the polynomi- 
als involved. For all cases we prove for the algebraic curve Ck(x, y) = 0, Ck(x, y) G 
Q[x, y], the following properties, which are similar to the ones given for k — 3. In 
each of the six cases there is a different value, say z G (—2, —1), which is given 
as a zero of a polynomial with rational coefficients computed through a suitable 
resultant. 

• The only multiple point is the origin, which is a double point. On it, the 
curve has two smooth branches tangent to two given lines with rational 
slopes; 

• the polynomial Ck(0, y) has the as a double root and k — 2 simple nonzero 
roots; 

• for each x G (z , 0) the polynomial y i — > Ck(x, y) has exactly k simple real 
roots; 

• the polynomial y i — > Ck(z , y) has k — 2 simple roots and a root of multi- 
plicity two; 

• for each x G [—2, z ) the polynomial y i — ► Ck(x, y) has exactly k — 2 simple 
real roots. 

□ 
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